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Abstract 

We investigate the shell structure of spherical nuclear bubbles in simple phe¬ 
nomenological shell model potentials. The bunching is produced by the energy 
gaps between single-particle states with a different number of radial nodes. The 
shell correction energies for doubly magic bubbles may be as large as -40 MeV and 
probably imply a very long lifetime against spontaneous fission, /t-stability occurs 
for ratios of the neutron number N to the proton number Z which differ markable 
from the /3-stability valley of ordinary compact nuclei. The a-decay probability is 
shown to be very small for proton rich bubbles with a moderately large outer radius. 
Metastable islands of nuclear bubbles are shown to exist for nucleon in the range 
A=450-3000. 

1 Introduction 

The idea that stable or metastable nuclear bubbles might exist is rather old. Already in 
1946 Wilson [Q has written about spherical bubble nuclei. J.A. Wheeler p|, in his famous 
notebook, mention the possibility that nuclei with different topology might exist. In 
1967, P.J. Siemens and H.A. Bethe [§] studied the energy of strongly deformed spheroidal 
nuclei and of spherical bubble nuclei using a liquid drop model (LDM). They looked at 
nuclei in the region of the then expected superheavy nuclei but did not investigate the 
stability of the spherical bubbles with respect to deformations. F. Wong was the first to 
perform detailed studies: Within the pure LDM he showed that for a fissibility parameter 
Xo > 2 spherical nuclear bubbles have a lower energy than compact spherical nuclei 
but are not stable against fission. He then investigated the shell correction energy for 
toroidal and spherical bubble nuclei on the basis of a harmonic oscillator potential using 
the Strutinsky’s method. In the study of the shell correction energy, he restricted his 
attention to known /3-stable nuclei and found spherical bubbles with a very small inner 
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radius as lowest energy configurations for certain doubly- magic known nuclei 

More recently, L. Moretto et al. [Q investigated the effect of an ideal dilute classical 
gas of nucleons at finite temperature inside the bubble and claimed that the pressure 
exerted by this gas on the layer of nuclear matter could stabilize spherical bubbles against 
deformation. As the authors pointed out, this mechanism may also be of relevance for 
the existence of ’’mesoscopic bubbles”. 

In what follows we investigate the effects of the shell structure for spherical bubbles in 
a large range of neutron (N) and proton [Z) numbers considerably extending them beyond 
the known nuclear species. The present paper is an extension of our study presented in 
Ref. 0. We use phenomenological shell model potentials and apply the Strutinsky 
method to evaluate the shell correction energy 6Esheii- The motivation of this work is to 
hnd the magic neutron and proton numbers and to obtain an estimate of the magnitude of 
the shell energy. This simple approach is expected to be sufficiently realistic to establish 
the importance of shell effects in nuclear bubbles and to lead to a hrst assessment of the 
stability of such systems. 

Surely, a reliable investigation of the existence of nuclear bubbles with an appreciable 
life-time versus hssion and other decay channels requires careful and technically difficult 
Hartree-Fock (HF) or Hartree-Bogoliubov (HB) calculations. Indications that bubble 
solutions may exist have indeed already been found recently in HB-calculations with 
Gogny interactions |0] and some time ago in HF-calculations based on the Skyrme 
interaction. 


2 Theory 

As in the original shell model of M. Goppert-Mayer and J. H. D. Jensen [^, we investigate 
the shell structure in the two limiting cases of a square well with inhnite walls, on the one 
hand, and an oscillator potential, on the other. 

In order to incorporate the case that the interior of the bubble is not entirely empty 
but hlled by an internal halo of nucleons, we also investigate the effect of replacing the 
above mentioned single particle potentials by a constant potential of (small) hnite depth 
V 2 for radial coordinates smaller than the inner radius {R 2 ) of the bubble. 

The spin-orbit potential for a bubble nucleus is smaller than for an ordinary compact 
nucleus because the contributions arising from the inner and outer surface have a different 
sign. Therefore, we treated the spin-orbit potential in perturbation theory. 
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2.1 Eigenfunctions and eigenenergies of the spin-independent 
potentials 

In the coupled representation which will turn out to be useful for the perturbative treat¬ 
ment of the spin-orbit potential, the eigenstates of the unperturbed Hamiltonian 

have the form 

yijm-= ^ {lrni^ms\l^jrn)Yirnt{0,(j))xrn, (2.3) 

Here, V(r) is a phenomenological spherically symmetric shell model potential. The quan¬ 
tities appearing in (2.3) have an obvious meaning. The radial wave functions (pni{f') are 
normalized eigenfunctions of the radial Schrodinger equation 


2_a 

2M r dr 


r2 j 


^nl + (y{r)- enl)(Pnl{r) = 0, 


which have to satisfy the boundary conditions 


(2.4) 


= 0 for all I 7 ^ 0; lirn(rv3„o(r)) = 0 , (2.5.1) 

r —^0 

LPni{oo) = 0 for all I . (2.5.2) 

The eigenvalues Sni and eigenfunctions depend on the orbital angular momentum 

I and on the number of radial nodes. Not counting the zeros at r = 0, oo, the number of 
radial nodes is given by (n — 1) where n = 1, 2,... . 
i) In the case of the inhnite square-well 


V{r) 


—Vq for R 2 < r < Ri 
-l-cxo otherwise 


( 2 . 6 ) 


the boundary conditions (2.5) are replaced by 


^nl{R2) = ^nl{Rl) = 0 . (2.7) 

A hnite well depth Vq (= 50 MeV) is introduced for later convenience (see Eq. (2.11)). 
It has no influence on the shell correction energy. The eigenstates (pnii'r) are linear com¬ 
binations of a spherical Bessel function ji{ar) and a spherical Neumann function yi{ar) 
(for ’’spherical Bessel function of the hrst and second kind”, see Ref. [p!0[| ) 

^ni{r) = Nni[ji{ar) + b^iyiiar)]. (2.8) 
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The parameter a is defined by 


a = 


l2M{eni + Vo) 


(2,9) 


and the normalization constants Nni are to be determined by the condition 

J drr^\ipni{r)\‘^ = 1 . 


( 2 . 10 ) 


The two parameters Eni and bni are obtained by satisfying numerically the boundary 
conditions (2.7). 


ii) In order to describe a finite density in 

the interior ir < R 2 ) 

1 of the bubble, we consider 

the slightly more complicated potential 




r -v^2 

for 

0 < r < i ?2 


f"(r)= -Id) 

for 

R 2 < r < Ri 

, (2.11) 

[ +00 

for 

r > Ri 



where the depth parameter V 2 (> 0) in the interior region is substantially smaller than 
Vo (> 0). As a rough estimate one may choose the ratio V 2 /V 0 to be equal to the ratio 
Pn{p)/ Pn{p) of fVe average density of neutrons (protons) in the region r < R 2 to the 
uniform LD density Pn{p) of neutrons (protons) in the region R 2 < r < Ri. The radial 

functions ipniij) and the eigenenergies Eni which correspond to the infinite square well 
with internal step (2.11) are given in Appendix Al. 

iii) As an opposite extreme to the one of the infinite square well we consider a harmonic 
oscillator with a finite internal step: 


V{r-) 


—V2 for 0 < r < i?2 (region II) 

—Vo + ^^^^(r — R)^ for r > i?2 (region I) 


( 2 . 12 ) 


where 

R '■= -^ 2 ) • (2.13) 

In Appendix A2 we exhibit the mathematical form of the radial eigenfunctions p>ni{f') for 
the potential (2.12) and discuss the numerical methods which may serve to obtain them. 

The WKB method is known to yield the eigenvalues of simple nuclear potential 
models rather well we have therefore applied it in our case. 

The upper (b) and lower (a) turning points for a particle of energy e and orbital 
angular momentum / in the potential (2.12) are given as the real solutions of the equation 


2Ma2 


V{a) 


nH{i + 1 ) 

2M62 


+ V{h ). 


(2.14) 
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As the analytic roots of this quartic equations are complicated, it is more convenient to 
solve (2.14) numerically. The WKB-eigenvalues Sni are dehned as the solutions of the 
equation 0 _ 

rb 


dr^e-V{r) 


r/(/ +1) 

2 Mr 2 


= {n + - ]h-K. 


(2.15) 


Here again an analytical integration is possible but impracticable whereas the numerical 
integration can be easily performed. 

So far, we neglected the average Coulomb potential Vcb(?’) acting on the protons 


Vcb(r) = eo / . ( 2 . 16 ) 

If one replaces the density distribution (r) of the protons by a uniform distribution in 
the layer between the inner and outer LD-radius 

p(p)(f) = 0o(r - i?2)0o(i?i - r)pp , (2.17) 


o 3Zeo 

47 r(fi? - fit) ’ 

the Coulomb potential Vchi'f') has the form 


Ccb(r) 


= f eo p, {e,{R2 - + do(r - Ri) 

+ do(r-i?2)do(i?i-r)(^-f-4)} . 


r 


(2.18) 


(2.19) 


If one wishes to take this potential into account, the potential V{r) must be replaced by 
\y{r) + Vcb(?")] whenever it acts on protons. 

We draw attention to the well-known result that a proton in the interior of the bubble 
(r < R 2 ) feels a constant Coulomb potential. This means that only the repulsive Coulomb 
force between protons in the interior region disfavours the formation of an internal proton 
distribution as compared to the formation of an internal neutron distribution. 

At this point, we would like to make the following remark: The shell correction 
energy is invariant with respect to a constant shift of all the eigenvalues of a given 

sort of particles. The shifts can be different for neutrons and protons. This is the reason 
why we may identify the nuclear part V (r) for n and p as long as we neglect the Coulomb 
potential and assume that the other parameters of the potential (radius parameter i?i, 2 , 
frequency a;) are the same for n and p. This simple assumption turned out to be em¬ 
pirically very successful for the Nilsson potential. Therefore, we adopt it too, realizing 
at the same time that it may be insufficiently accurate for cases of very exotic nuclear 
compositions. 

The effect of a hnite potential V 2 7 ^ 0 in the interior region r < R 2 is qualitatively 
the same for the two potentials (2.11) and (2.12). Therefore, we studied the effect of a 
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finite value of V 2 only for the infinite square well. In all calculations with the harmonic 
oscillator potential (2.12) the parameter V 2 was chosen to be 0 : 

V(r) = —Vo + (r — Ry for r>0. (2.20) 

Furthermore, we determined the eigenvalues Eni corresponding to the potential (2.20) from 
the WKB approximation and the radial s.p. densities 

Pni(r) = (fhir) 

from the Thomas-Fermi model. The details of this procedure are presented in Appendix 
A3. 

We still have to state how the ’’shape parameters” of the shell model potentials are 
to be chosen: Since we use the Strutinsky method representing the total binding energy of 
the nucleus as a sum of the liquid drop (LD) energy i^Lo and the shell correction energy 
(5i?sheih fhe shape parameters of the shell model potential must be chosen in accordance 
with the ones of the liquid drop. When calculating the energy as a function of the nuclear 
shape, the multipole moments (of low order) calculated with the shell-model wavefunction 
must agree with the ones calculated with the LD-density. For the case of spherical bubble 
nuclei, the MS deviation from the central sphere of radius R must be the same when 
calculated with the shell-model wave function and with the LD-density 

((r-i?)^)sM = ((r-^)^)LD, (2.21) 

where 

R:=^(R, + R2). (2.2F) 

We restrict our consideration to bubble nuclei with a negligible nuclear density inside of 
the smaller radius R 2 of the LD-distribution. This means that we only consider the cases 
of the inhnite square well and of the shifted oscillator without an internal step (see Eq. 
( 2 . 20 )). 

For the shifted oscillator, the condition (2.21) may serve to hx the frequency u. For 
the square well potential it is physically reasonable to choose the potential walls at a 
slightly larger spacing than the difference of the LD-radii 

Ri = Ri + AR , (2.22) 

R2 = R2 + AR . (2.22') 

The condition (2.21) may then serve to determine the parameter AR. One could modify 
these constraints slightly. So the shift (2.22) of the inner and outer radius of the potential 
with respect to the ones of the liquid drop need not be equal and the reference radius R 
(see (2.21’)) might be chosen differently. We believe that the simple choices we made are 
reasonable given the fact that the Strutinsky method itself is of limited accuracy. 
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The MS-radius calculated within the LDM is given by 


((r - ^ 


- i?2^) + 2{-Ri^R2 + Ri^R2^ - Ri^R2^ + RiR2^) 


{Rl^ - R2 


(2.23) 


The MS-radius ((r —i?)^)sM in the shell model is approximately calculated in the Thomas- 
Fermi approximation (see Appendix A3): 


{(r - fl)2)sM d?p (r - Bf [e,[el - - l/(r)] + 

+ Do [4 - i - r(r) - Vtar)]} . 


(2.24) 


The Fermi energies for the neutrons and protons are determined by the number N 

of neutrons and Z of protons 


Z = It / dV / rf3p9o [e^F - 4 - v{r) - Vct,(r) 

= /;J 14 - V{t) - l/cb(i')]»'2 . 

All integrations but one can be trivially done and we obtain 




tsM « + Jr) . 


(2.25.1) 


(2.25.2) 


(2.26.1) 


Jn := / drr‘^{r - Rf[e'},-V{r)]'^ , 


(2.26.2) 


Jp\= / (irr^(r - i?)^ - ■F(r) - V'cb(r)]' 


(2.26.3) 


The limits of the integrals are the (lower and upper) classical turning points of a particle 
with energy ep in the potential V (r) 

V(r^) = V{r^) = e^p (2.27.1) 

f"(rf) + Vck{4) = V{r{) + Vcb{r^,) = 4 . (2.27.2) 

We neglect the Coulomb potential in the numerical calculations. 

2.2 Spin-orbit potential 

In the conventionally used Skyrme interactions, the effective spin-orbit interaction be¬ 
tween two nucleons has the form 


uso := i Wo{di + a 2 )k x 5{ri - f 2 )k , 


(2.28) 
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where Wn ~ 120 MeV fm^ and 




V 2 ). 


This interaction leads to the following one-body spin-orbit potential ||I2 


(2.28') 


^so = fh(r)-(vxcf) , 
% 

fy(r) := —— VP(n • 


(2.29) 


(2.29') 


Here p{f) is the sum of the density distribution of neutrons (pn(fO) protons (pp(T)). 
Recently, Reinhard and Flocard generalized the ansatz (2.28) for the spin-orbit in¬ 
teraction in such a way that the vector W{r) (2.29’) takes the more general form 


fRn(p)(r) = h VP(^) + K ypn{j,){r) 


(2.30) 


for neutrons and protons, resp. Here, 64 , 64 are parameters in the notation of Ref. 
which have the values 64 = 64 ,= 62.130 MeVfm® for the Skyrme interaction Skll. This 
choice corresponds roughly to the value of the parameter Wo ~ 120 MeVfm®, if the 
numbers of neutrons (N) and protons (Z) are equal. For unequal numbers of n and p, 
the form (2.30) yields are more flexible form of the spin-orbit potential. With the choice 
64 = 64 = 62.130 MeVfm® one reproduces the empirical values of the isomer shifts across 
magic shell closures [^. 

For spherically symmetric density distributions the spin-orbit potential Vso takes 
the form 


fyn(p) _ 
•^SO — 




dr 


dr 


I ■ s, 


(2.31) 


where s is the dimesionless spin-operator and I is given by 


I := fx -V • 


(2.31') 


For symmetric spherical nuclei 


pn{r) =Pp{r) = -p{r) 


one obtains the still simpler form 

_ T/(p) _ 


*^80 — •'so — •'SO — ^ (t»4 -h ^0^) L ■ S , 


(2.32) 


(2.33) 


which is equivalent to the conventional isospin independent form of the spin-orbit potential 
of a spherical nucleus (see Eqn. (5.103) in Ref. [11]) 
















So far, our calculation of the spin-orbit coupling is based on the isoscalar form (2.33’) of 
the spin-orbit potential. Furthermore, we used phenomenological forms for the densities 
p„(p)(r) rather than the ’’selfconsistent” densities which could be obtained from the single¬ 
particle functions. The following density distributions were used 


PnUr) = I --( 2 . 34 . 1 ) 

Pn(p)(r) = p:<^> { (l + e"^) (l + e^) }"‘ ( 2 . 34 . 2 ) 

Pn(p)(r) = - fl 2 )«„(fli - r)(l - K(r - Rf) . ( 2 . 34 . 3 ) 

where R = + i? 2 )- The parameters are determined by number conservation 


N(Z) = J d^rpn(p){r) . 


(2.35) 


The results obtained from the different forms of Pn{p){f') are very similar. The 
parameter k in (2.34.3) is chosen in such a way that the density at r = Ri( 2 ) is equal 

Po^"V2 : 


2 

{Ri - R2Y ■ 


(2.36) 


The parabolic choice (2.34.3) of the density yields a particularly simple form of the spin- 
orbit potential: 


14 


n{p) 


SO 


12fFop;)(p) (r-R\f - 

{R, - R^y ) 


(2.37) 


Using the unperturbed eigenstates in the coupled representation (2.2), we obtain the shift 
^Snij of unperturbed eigenvalue Sni due to the spin-orbit potential 


= -ITo ■ 47r 


/•“ dp{r) o 

/o 


j(j + l)-/(/ + l)-^ 


(2.38) 


^nlj ^nl • (2.39) 

In Eq. (2.38), the function p„; (r) is the unperturbed single particle density 

°Pnl= J = <4nz(^)<4nz(r) . (2.40) 

The total angular momentum j in (2.38) can assume one of the two values j = |/ ± ||. 
The approximate calculation of p^i (r) is described in Appendix A3. 
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3 Results 


3.1 Bubble nuclei in the LDM 


As we have stated already in the introduction, C.Y. Wong g] investigated the stability 
of bubble nuclei within the liquid drop (LD) model. Also the role of the shell effects has 
been hrst discussed in Ref. [^. The result was that while for large enough mass (A) 
and charge numbers (Z) spherical bubble conhgurations do have a lower energy than the 
compact spherical ones, the energy can still be lowered by deforming the bubble. This 
means that no barrier protects the spherical bubble nuclei from undergoing hssion. 

It is useful to represent these results in a simple form. We denote the inner and 
outer radius of the spherical bubble nucleus by R 2 and Ri, resp., and the radius of the 
compact spherical nucleus of the same mass number by Rq. We retain the conventional 
basic assumptions of the LDM that the charge is uniformly distributed over the nuclear 
matter and that the nuclear density is constant inside the matter distribution dropping 
abruptly to zero at the surface. 

The difference AZld between the energy of the bubble conhguration with radii Ri ^2 
and the energy of a compact spherical nucleus of radius Rq is given by 

AZlo = 4:Tia{Rl + RI — i?o) + EchiRi, R 2 ) — T^cb(Ro), (3.1) 


where a is the surface tension parameter. 

The Coulomb energy Zcb(Ri, R 2 ) of the spherical bubble nucleus is dehned as a func¬ 
tion of the proton density pp{r) and the Coulomb potential Vc'fe(n) by 

Ec^{Ri, R 2 ) = ^J (fr Vch{r)pp{r). (3.2) 

The LD-density Pp{r) is represented by the simple formula ( 6^0 = Heaviside function) 

Pp{r) = 9 q{Ri - r)6'o(r - R 2 ) Pp (3.2') 

and the corresponding Coulomb potential is given by Eq. (2.19). The Coulomb energy of 
the LD-bubble is then obtained in the form 


Ecb(Ri,R2) = ^^ 


^ p5 3 p5 _ p2 p3 

III T cr-^2 -^1-^2 

5 5 


o 3Zeo 

The Coulomb energy Ech{Ro) of the homogeneously charged sphere is 

3Z2e^ 


Ech{Ro) — 


5Rn 


(3.3) 
(3.3') 

(3.4) 


The condition of volume conservation yields the relation 


rp3 _ p3 p3 

JTq — K2 


(3.5) 
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It is easily seen that the energy difference Ai^LO in units of the surface energy of the 
spherical compact nucleus dvrcri^o can be written as 




fl(f2) + |f2 - |>’2fl(f2) - 1 


(3.6) 


Vi and V 2 are the radii of the outer and inner surface in units of Rq, resp. They are related 
by: 

p / p3 \ 1/3 

-^1 / 1 , -^2 \ 

Ro 

Xq is the conventional fissility parameter defined by 


-I=^=(I+= (I+ 


(3.7) 


A'„ = 


-E'cb(.Ro) 


3Z^el 


(ZVA) 


2 ■ AnaRl AOnaR^ {ZyA)^nt ' 

The dimensionless parameter {Z‘^ j defined by the last Eq. (3.8) 


(3.8) 


(ZVA)„it = 


dOvrarQ 

3e§ 


(3.80 


depends on the value of the surface tension constant a and the radius parameter tq. For 
ro = I, 2 fm and cr = 1.026 MeV fm“^ one obtains {Z‘^ j = 51.57. We draw attention 
to the fact that the difference (in units of the surface energy Es = dvrcri^o) ^ between the 
energy of the bubble configuration and the energy of the compact sphere depends only 
on the fissility parameter Xq. 

The condition of stationarity reads 


dV2 


2 u| 

Vl{v2) 


+ 2 v 2 + 2Xo 


5 

2 


vivl{v2) + yUa 


5ui 

Vl{v2) 


0 . 


(3.9) 


where vi{v 2 ) is defined be Eq. (3.7). 

Beside the trivial solution U 2 = 0 (compact sphere), the equation (3.9) has 2 addi¬ 
tional real solutions^ if and only if Xq > 2.02. We can write Eq. (3.9) in the form 



V2+Vl{v2) 

V2 3^2 (U2 - Vi{v2) 


(3.10) 


The inverse U 2 (Xo) of (3.10) describes the values of V 2 which correspond to stationary 
values of the energy S. The function is shown in Fig. I. From this figure, one sees in an 
illustrative way that for given Xq > 2.02 there are 2 real positive solutions V 2 {Xq). As is 
to be seen in Fig. 2, the smaller one (u^(Xo)) corresponds to a maximum and the larger 
one (u^(Xo)) to a minimum of the energy S. 

iReal solutions with still larger V 2 may exist but they are physically meaningless. 
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For small values of V 2 (^2 1)) the energy change S by bubble formation can be 

written approximately as 

+ 5Xo)vl + 3Xovl + 0{vl). (3.11) 

This shows that the compact spherical LD (^2 = 0) is always ’’locally” stable with respect 
to bubble formation independently of the value of Xq. In other words: If there is at all 
a bubble solution, then it is separated from the compact sphere by a barrier. 

It is a curious result that the inner and outer radius of the bubble nucleus (measured 
in units of i?o = depend only on the hssility parameter Xq. If the parameter 

{Z^/A)crit (see (3.8’)) is independent on (X — Z), a given value of the hssility parameter 
Xq implies a given value of Z"^/A. 

For (Z^/y4)crit=51.57, bubble solutions within the LDM are only obtained for 


Z^ 

— > 2.02 X 
A 



(3.12) 


Thus, for given nucleon number A, the charge must be larger than a minimal value for 
LD-bubble solutions to exist. From (3.12) we obtain Z > 250 for A = 600 and Z > 353 
for A = 1200. The minimal proton density, necessary for LD-bubble solutions to exist, is 
thus the larger, the smaller the nucleon number A. 

Let us restate clearly://we speak about the existence of a bubble in the LD theory, we 
only mean that a minimum of the LD-energy with respect to (volume-conserving) changes 
of the bubble radii exists. As we shall show later, and as it was shown by the author of 
Ref. 0, these spherical solutions are unstable against deformation of the bubble. 

Nevertheless, it is physically meaningful to study carefully the conditions for the 
existence of bubble solutions in the LDM, because whenever these solutions exist, the 
shell structure effects have a better chance to make the bubbles stable. 

The Figs. 3-5 are designed to give a global view on the LD-properties of bubble 
nuclei. There results were obtained with the parameters of the liquid drop model taken 

0 . 


from Ref. 


In Fig. 3, lines of constant hssility parameter 


Ech{Ri, R2) 

87ra(Rf + Rl) 


(3.13) 


are shown in the (X, Z)-plane. Each point corresponds to an existing bubble solution. 
We draw the reader’s attention to the fact that the hssility parameter X dehned in Eq. 
(3.13) dihers from the hssility parameter Xq, which refers to a compact spherical nucleus 
of radius Rq. The two hssility parameters are related by 


X = X, 


n? + 


3,,5 
2 ^2 


2 ^2 


t’? + v'i 


= A'„ 


(1 + - If) 

(1 - - / + / 2 / 3 ) 


(3.14) 
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where the parameter / is dehned as the fraction of the empty volume to the total volume 
of the bubble 

/=(!)'. (3.15) 

Contrary to Xq, the fissility parameter X depends on the two reduced radii of the bubble. 
It is a measure of the hssility of the spherical LD bubble. The larger the value of X the 
more rapid is the descent of the energy as a function of the deformation. 

Each point on the contour lines shown in Fig. 3 corresponds to a bubble solution. 
The parallel straight lines represent isobars. As one moves on a given isobar towards 
nuclei with a larger proton density, the fissility parameter in general decreases, i.e. the 
LD-bubbles become less unstable versus hssion. A peculiar feature is that the line of 
hssility X=1.25 is cut twice by the isobaric lines. This means that for given A there 
may exist two bubble solutions corresponding to the same hssility X but diherent charge 
numbers and diherent bubble radii. Of course, the 2 solutions belong to diherent binding 
energies in general. 

Lines of equal binding energy gain AEld (see (3.1)) in the {N, Z)-plane are shown 
in Fig. 4. The parallel straight lines indicate again isobars. On the proton rich region of 
the hgure for given A, the binding energy gains due to bubble formation are the larger, 
the larger the charge Z. This can be easily understood by the fact that the decrease of the 
repulsive Coulomb energy due to bubble formation is larger for larger charges. Certain 
equi-energy lines are cut twice by given isobars. This means that there are two isobaric 
nuclei {Ni, Zi) and {N 2 , Z 2 ) with bubble solutions corresponding to the same energy gain, 
but in general diherent bubble radii. Following the isobaric lines in the direction of larger 
binding energy, one can infer the type of /5-decays to be expected. Curiously, on the p- 
rich side, /5+-decays would increase the binding energy even further. These considerations 
must of course be taken with a grain of salt, because the bubbles are in general not stable 
with regard to hssion. Stability can only be produced by shell ehects which complicate 
the topology of the equi-energy lines. 

In Fig. 5 we display lines of equal hole parameter / in the {N, Z)-plane. Large / 
means that the interior hole represents a large part of the total bubble volume. It is seen 
that in general, as one increases the charge Z along an isobaric line, the hole fraction 
increases. Again there are pairs of isobaric nuclei {Ni,Zi), {N 2 , Z 2 ) with bubbles of the 
same shape at diherent neutron and proton numbers. The energies of such pairs are 
usually diherent. 

In Fig. 6, we display lines of equal total binding energy per particle for LD-bubbles 
in the {N, Z)-plane. It is remarkably diherent from Fig. 4 where we showed lines of 
constant total binding energy gain due to bubble formation. As a general trend, the 
total LD-binding energy per particle tends to decrease as the proton number increases 
along an isobaric line. The largest binding energies per particle are seen to occur in 
the region 1000<N<1800 and 300<Z<400. It is also interesting to see an island of very 
small (absolute) binding energy per nucleon in the region of 400<N<600 and 380<Z<480. 
Again we observe pairs of isobars showing the same binding energy per particle. 

Let us emphasize that the simple features of the Figs. 3-6 are due to the simplicity 
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of the LDM. As a macroscopic model it predicts a smooth dependence of all the physical 
properties on N and Z. The shell effects, which will be shown to prevent certain magic 
bubble nuclei from undergoing hssion, will complicate the picture considerably. The shell 
structure will be seen to modify the smooth trends of the LDM profoundly, but only in 
certain localized regions of the {N, Z)-plane. 

So far, we only considered spherical nuclei. A crucial question is whether the spher¬ 
ical bubbles are stable with respect to deformations of the inner (S' 2 ) and outer ( 51 ) 
surface. Within the LD approximation, we have to study the energy difference 


Ah;LD(^2, ^ 1 ) := ^Ld(^2, ^ 1 ) - ^LD(i?2, ^l) , 


(3.16) 


where T^ld(*S' 2 , S'!) is the energy of a deformed LD bubble with the inner and outer surface 
S '2 and Si, resp., and Ei^d{R 2 ,Ri) is the LD-energy of the spherical bubble with radii 
R 2 , Ri- 

As the volume terms cancel again, we have 




= (T [S'2 -|- S'! — AttRI — 
+ D J PpTi)PpT2) ^3 ^3^ 

[2 •> \rr—r2\ ^ 


AttRI] 

2 — Eci,{R 2, Ri) , 


(3.17) 


where S' 2 (i) denote the magnitude of the deformed inner (outer) surface and Pp{r) is the 
uniform charge distribution within the volume enclosed by S '2 and S'!. An investigation 
of AA'ld shows [^ that the spherical LD bubbles are not minima but saddle points in 
a multi-dimensional potential surface. In particular, the energy of a spherical LD bubble 
is lowered by quadrupole deformations. 

In Fig. 21 we show the energy AF^ld as a function of the quadrupole deformation 
of the outer surface. The energy is seen to decrease for prolate as well as for oblate 
deformation of the bubble. It is also seen from this hgure that constraints on the shape of 
the inner surface (for given deformation of the outer one) increase the energy substantially. 

As is well-known, ordinary compact spherical liquid drops are unstable with respect 
to hssion for hssility parameters Xq > 1. This means that the formation of a bubble 
nucleus by a fusion of 2 nuclei, is complicated by the fact that only a tiny part of the 
reaction cross section is expected to contribute to the formation of a bubble. 


3.2 The level spectrum 

As explained in Sect. 2, we calculated the level spectrum for the shifted inhnite square 
well and the shifted oscillator treating the /s-coupling in perturbation theory. We also 
studied the effects of a hnite negative constant potential in the interior region for the case 
of the square well. This corresponds to the case of a bubble nucleus with a reduced but 
nonvanishing internal density. 

In Fig. 7, we show the level spectrum for the shifted inhnite square well as a func¬ 
tion of the hole fraction /. A nice feature of the shifted inhnite square well is that its 
eigenvalues scale with We have used the unit 100/A^/^ MeV. Consequently, from 

the spectrum given in Fig. 7 one can infer the eigenvalues for arbitrary nucleon number A. 
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We use the conventional spectroscopic notation {nlj)\ n = 1, 2,... counts the eigenvalues 
for given /,j in rising order. Apart from s,p,d (/ = 0,1,2) the angular momentum / is 
represented by the letters in alphabetical order {f,g,h,... ~ / = 3,4, 5 ...). 

The following physical features of the spectrum should be noted: for n = 1, the single 
particle energies Cnij decrease as a function of /, for n > 2, they increase as a function 
of / with a steepness which becomes drastically larger with n. The dependence on n is 
produced by the requirement of a vanishing wavefunction at the limits of the well. As / 
increases, the diameter d = Ri — R 2 diminishes and, consequently, the energy difference 
between successive radial modes increases. 

The strong dependence of the difference between successive eigenvalues of the same 
orbital angular momentum I but different radial modes can be understood in the following 
simple approximation: If we replace the centrifugal term radial Schrodinger 

equation by its value at the center R = of the bubble layer, the eigenfunctions 

Unii'f') become trigonometric functions which are either ’’even” or ’’odd” with respect to 
the reflection [r — R) ^ —{r — R). In this approximation the eigenenergies have the form 
(n' = 0,1,...) 


h^/(/ + 1) f (2n' + 1)^ for even states 

2MR? I (2n'+ 2)^ for odd states 


One can easily check that this relation explains qualitatively the behaviour of the eigen¬ 
values as a function of the hole fraction /. We note that 


d^ = Rl 


R^ = 


(1 _ /l/3)2 

(1 - /) 2/3 


4 (l-/p/3 

and that the spectroscopic label n is related to n' by 


n = 

1 ~ 

n' 

= 0, 

even state 

n = 

2 ~ 

n' 

= 0, 

odd state 

n = 

3 ~ 

n' 

= 1, 

even state 


a.s.o. 


(3.16.1) 

(3.16.2) 


(3.17) 


The gentle dependence of the energies with n = 1 as a function of / is due to the fact 
that in this case the increase of the second term in (3.15) is largely compensated by the 
decrease of the hrst term in (3.15). 

Let us stress that the approximation (3.15) is crude and cannot replace the deter¬ 
mination of the true unperturbed eigenvalues by a numerical solution of the Eqs. (2.7). 
In Fig. 7 we plotted the energy levels including the spin-orbit interaction 

Cl/ (^nlj ^nl T (t/’nZjr'm I LsO I• (3.18) 

An unusual feature is that the levels with j = / — ^ lie below the levels with j = I + 
Furthermore, the spin-orbit splitting, i.e. the energy difference is much 
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smaller than in the case of ordinary nuclei with the same 1. The reason for both these 
phenomena is difference of the density distributions: Since the spin-orbit term depends on 
the gradient of the density (see Eq. (2.29’),(2.30)), its main contribution originates from 
the nuclear surface. The outer surface contributes with the same sign, but the inner with 
the opposite sign as compared to ordinary nuclei. For spherical bubble nuclei the spin- 
orbit potential can be written in the form (2.31). The factor 1/r in the expression weighs 
the inner surface more than the outer one so that the sign change becomes understandable. 
The inhnite square well model is, of course, not realistic especially for bubble nuclei with 
a small value of /. In the case of bubbles with a small hole fraction, we cannot expect 
the hole to be empty. It will rather be hlled with a reduced nucleon density and thus the 
change at the inner surface from the density in the layer to the reduced density in the 
hole region will be much smoother than in the case of the inhnite square well. In this 
case, the contribution of the outer surface may become more important than the one of 
the inner surface in spite of the weight factor 1/r and thus the usual order of spin-orbit 
splitting ensues. 

In the Fig. 8 we show the level scheme as a function of the parameter / for the 
shifted oscillator. The oscillator is more suitable than the inhnite square well (i) for small 
/-values, where the inner region is expected to contain a reduced density of nucleons, 
and (ii) for bubbles with a thin layer, i.e. with a diameter d which becomes smaller than 
about 4 fm. As one can infer from (3.16.1) this happens for large values of / and at the 
same time not too large values of A. (Example: /=0.3 (0.4) yields d/RQ=0.37 (0.31) and 
consequently d=3.9 fm (3.3 fm) for A=700.) Contrary to the case of the inhnite square 
well, there are also cases of ’’normal” order of the spin-orbit splitting (see for instance the 
pairs (^2iY,2iY^ and (i/^y; 1^^)) harmonic oscillator potential. 

The following general trends can be stated: 

(i) For increasing hole fraction /, the contributions of the inner and outer surface 
region to the total spin-orbit term tend to become more equal. Since the two 
regions contribute with opposite sign, the spin-orbit splitting tends to decrease with 
increasing /. This is seen in Fig. 7 and 8, i.e. for both the potentials. 

(ii) The centrifugal potential has the tendency of shifting the single particle density to 
larger values of /. For given hole fraction /, this ehect increases with I and for given 
/, the ehect decreases with /. This shift of the s.p. density towards larger values of 
r is expected to increase the contribution of the outer surface region as compared 
to the inner one as far as the spin-orbit term is concerned. 

In fact one can observe in Fig. 8 that for large enough I and small / the ’’normal” 
spin-orbit splitting prevails. The ehect is more pronounced for the oscillator than for the 
square-well. 

In Fig. 9 and Fig. 10 we present the level scheme for the square well potential 
(2.11) with a hnite potential step (14 — Vq) = 60 MeV and 40 MeV, resp. The hnite 
internal potential step presents the possibility to study qualitatively the ehect of a hnite 
reduced density inside the bubble. In a classical description this was done in Ref. 0. An 
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internal density distribution changes, of course, the surface tension of the LD at the inner 
surface and, to be consistent, we should also add a LD energy describing the macroscopic 
smooth energy originating from the internal density. The results of Fig. 9 and 10 are to 
demonstrate the effect of a reduced internal density on the level scheme. For this purpose 
we compare the level scheme in Fig. 9 and 10 with the one in Fig. 7: 

(i) The density distribution extends from the inner surface down to r = 0 contrary to 
the case of the inhnite square well of Fig. 7. Consequently, the weight of the inner 
contribution to the spin-orbit splitting increases as compared to the inhnite square 
well. This is seen by the fact that the splitting in the unusual order tends to be 
larger in Fig. 9 and 10 as compared to Fig. 7. 

(ii) For energies far below the step, the results are very similar to the ones of the inhnite 
square well. 

(hi) For energies above the step, the changes are of course very large the general tendency 
being a compression of the spectrum at energies above the step. 

We would like to add a remark concerning the neutron (or proton) numbers which 
appear in the plots for the 3 values: 0.12, 0.24, 0.28 of the hole fraction parameter /, 
whenever the distance Siy+i — to the next higher s.p. energy exceeds 1.5 energy units. 
They represent the number of neutrons (or of protons) occupying s.p. states with energy 
Ek < Sy Only if the total energy of the bubble nucleus has a minimum at one of these 
/-values, we may expect that the numbers given represent magic numbers for bubble 
nuclei. Nevertheless, they are useful to give us a hint as to where magic shell closures 
might occur. Clearly, the precise values of the magic numbers depend on the type of shell 
model used. 

An interesting aspect can be observed concerning the dependence of the ’’magic 
numbers” on the hole fraction variable /. A general feature of all the level schemes as 
a function of / is the interplay of the steeply rising energies Cnij with n > 1 and the 
smoothly /-dependent energies euj. Magic numbers which appear in a region of many 
crossings depend sensitively on /, whereas magic numbers which occur between 1/j-levels 
away from bunches of strongly rising 21 j- and 3/j-levels depend only weakly on /. The 
strongly /-dependent magicity occurs preferentially for small values of /, and the weakly 
/-dependent magicity for large / values. 


3.3 The shell correction energy 

From the spectrum of s.p. levels we calculate the shell correction energy in the 
well-known way [13] (z/ = nlj] q = 1 : neutrons, q = 2 : protons) 


6E,^,niN,Z-,f) = E / ' ^9<iie)de - E 




<?=! 


(3.19) 
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(3.20) 


where gq{e) is the average level density: 



and fcorr is a 6th order correction polynomial of the variable = (e — • 


/corr — 1 + (2 ~ + (g ~ ^16 ~ 


(3.21) 


The position of the Fermi energy (A^) for the smooth energy spectrum is obtained by the 
integral: 

(or Z) = f gq{e) de . (3.22) 

J —00 

The ’’plateau condition” with respect the smearing width 7 was reasonably well fulhlled. 

In Fig. 11, we show the shell correction energy for the inhnite square well as a func¬ 
tion of the number of protons or neutrons. The hole fraction parameter is chosen to be 
/=0.24. The hrst substantial shell effect is seen to occur for N (or Z)=242. It corresponds 
to a negative shell energy of about —15 ■ 100/A^/^ MeV . For a doubly magic nucleus 

the shell energy may amount to (—25 to —30) ■ lOO/A^/^ MeV . The magic numbers and 
the magnitude of the shell energy depend on the hole fraction / (see Figs. 17-19). 

In Fig. 12, the dependence of the shell energy on the height of a constant potential 
in the internal region is studied: If only low-lying levels are occupied at a substantially 
lower energy than the step, the internal step has little influence. This is so for particle 
numbers < 200. For larger numbers of N or Z, the magic numbers and the magnitude of 
the shell energy differ noticeably from the case of the inhnite square well. In particular, 
the strong shell effect for 242 disappears completely whereas a new region of negative 
shell energy emerges at 330. 

For comparison we show the shell energy for a shifted harmonic oscillator corre¬ 
sponding to the same hole fraction / = 0.24 in Fig. 13. The frequency of the oscillator 
was chosen for the particle number A=500 in Fig. 13. Strong shell effects are seen to 
occur for several numbers of particles of a given sort. They are not the same numbers 
as for the square well (of the same /), as one would expect. Nevertheless, corresponding 
shell closures are not very far away from each other and the order of magnitude of the 
shell energy is the same. 

The dependence of the shell energy for given (iV, Z) on the hole fraction / is shown 
in the 3-dimensional plots of Fig. 14 (for an inhnite square well) and Figs. 15 and 16 
(for a harmonic oscillator corresponding to yl=500 and A=1000, resp.). It is seen that 
the magnitude of the shell energy increases with N (or Z). Each of the different islands 
of large negative shell energy corresponds to specihc limited intervals of the hole fraction 
/. These intervals are of different length A/ depending on the specihc magic number. 
On the average, the length A/ 0.05. The fact that magic neutron or proton numbers 
for nuclear bubbles are correlated with given hole fraction values / and the fact that we 
expect the n- and p-distribution in a given bubble nucleus to exhibit the same shape. 
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limit the possible double magic combinations. Doubly magic numbers for neutrons and 
protons can only occur in a given bubble nucleus if they belong to the same value of the 
hole fraction /. This is analogous to the shell effects of ordinary nuclei as a function of 
the quadrupole deformation. 

The correlation between particle number [N or Z) and /-parameter, which is visible 
in the Figs. 14-16, is shown once more in the form of undistorted 2-dimensional maps in 
Figs. 17-19. 

For ordinary exotic nuclei, especially for nuclei with a large excess of neutrons over 
protons, the distribution of neutrons and protons do not always coincide. In fact the 
radius of the n-distribution is probably systematically somewhat larger than the one of 
the p-distribution (”n-skin”). For bubble nuclei with a large n-excess we have to expect 
the same phenomenon. Thus, for a very n-rich bubble nucleus, we expect that the inner 
(outer) radius of the n-distribution is smaller (larger) than the inner (outer) radius of the 
proton distribution. This relaxes somewhat the condition of equal /-parameter for n and 
p. We have neglected this possibility in our present investigation. 

Let us now discuss the total energy: 

E = F'ld + ^T^sheii, (3.23) 

where F^ld is the liquid drop energy with a suitable choice of the reference energy. So far, 
we studied separately the dependence of F^ld and SEs^eii on the hole fraction / for given 
N and Z. For this purpose we used the LD energy of a compact spherical nucleus (radius 
Rq = as a reference energy (see Fig. 4) 

AEld := E^^{f- N, Z) - F;ld(/ = 0; N, Z) (3.24) 


because in this way the volume energy cancels due to the conventional saturation assump¬ 
tion. 


If one adds SEsheW to AT^ld one hnds bubble solutions in certain regions of the 
{N, Z)-plane. This is shown in Fig. 20. 

The crucial question is the stability of these spherical bubble conhgurations with 
respect to deformations. A careful investigation of the stability versus volume conserving 
deformations of the inner {S 2 ) outer (Fi) surfaces of the bubble nucleus requires the cal¬ 
culation of the surface- and Coulomb energy of an arbitrarily shaped bubble and — what 
is much more difficult to achieve — the calculation of the eigenvalues of a correspondingly 
deformed single particle potential. 

We have not yet performed such a calculation. It seems also a bit questionable to us 
whether such a big effort is meaningful in the framework of the Strutinsky method, given 
the fact that the Hartree-Fock theory (HF) and Hartree-Bogoliubov theory (HB) provide 
a quantitatively more reliable approach^. Instead, we applied a simple phenomenological 
method which had been introduced by W. Swiatecki in 1966. The ansatz is designed 


^The existence of stable bubble configurations for various mass- and charge numbers is being inves¬ 
tigated within the HB-theory on the basis of the Gogny-interaction by J. F. Berger and J. Decharge. 
Encouraging results were obtained by these authors (priv. com.). 
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so that the shell energy of the spherically symmetric nucleus is damped away as a function 
of the quadrupole deformation. The result is shown in Fig. 21 for a bubble nucleus of 
A=750 and Z=288 and a spherical bubble solution at /=0.26. The LD-energy at the 
spherical bubble shape is put equal to zero. The following features are of interest: 

(i) The LD-energy decreases for positive as well as for negative values of the quadrupole 
deformation (32 of the outer surface Si. 

(ii) Constraints on the shape of the surfaces (same deformations for the two surfaces or 
spherical inner surface S' 2 ), results in noticeably higher energy than minimization 
of the LD-energy for given {32 with respect to a limited number of shape degrees 
of Si and S' 2 . Even spurious flat minima on the oblate side are produced by these 
constraints. 


(iii) (5Esheii + -Eld exhibit a deep and relatively narrow valley with a minimum of about 
-30 MeV at (32 = 0. The shell effect disappears at about /3 ~ ± 0.2 using Swiatecki’s 
phenomenological form of (fEgheii- 


(iv) The barrier to hssion is seen to be very thick. The spontaneous hssion lifetime can 

safely be assumed to be practically inhnite. 

We can be sure that shell energies 5Esheii < -(20 to 30) MeV protect the spherical 
bubble nuclei from undergoing hssion, i.e. lead to a practically inhnite hssion lifetime. 

An important question is the dependence of our results on the parameters of the 
model: In the Strutinsky method, the surface tension cr and the radius parameter tq 
appear as phenomenological quantities. Apart from these two LD-parameters, which 
enter in a signihcant way (see Eq. (3.8’)), the form of the phenomenological shell model 
and the parameters therein inhuence the results. We note that the dependence of the 
surface tension a and the radius parameter tq on ((iV — Z) /A) ^ is of great importance for 
the stability regions of nuclear bubbles. Whereas we have some good empirical evidence 
on the isospin dependence of ro [0 , very little is known about the isospin dependence of 
a. 

In order to test how the existence of bubble solutions depends on the parameters 
of the nuclear models we also performed a calculation based on Skyrme’s energy density 
functional E[p\ for a spherical nucleus [|^. We have assumed that the kinetic energy 
density of protons Tp and neutrons are functions of the corresponding nucleon densities 


7 rV 3 .35/3 

Tq{r) = --- Pqi'f) were q = p,n, (3.28) 

i.e. we have neglected the higher order correction which appears in the extended Thomas 
Fermi approximation. We represented the densities Pni'^) and Pp{r) of n and p by Fermi 
functions localized around R = ^{Ri + R2) 


Pq{r) =pq 


-I I r—R\ H ^-^2 

\l + exp—^ 1-f exp—^ 


(3.29) 
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and determined the parameters Rl,R 2 , and the snrface thickness by variation of E[p]. 
Choosing the parameter set SkP proposed in Ref. [0 and considering a nnclens of 


R=470 and Z=184 we fonnd the resnlts shown in Fig. 22-25. To onr snrprise we obtained 
a beantifnl bnbble solntion for a nnclens with twice the nncleon and proton nnmber of 
92 ^Ui 43 ! This resnlt is not obtained for the more conventional Skyrme parametrization 
Skill. The specihcity of the parameter set SkP is that the effective mass of the nncleon 
is eqnal to its free mass. This £t leads to a more realistic level density at low excitation 
energy than other Skyrme parametrizations. Here we are not pleading for or against the 
interaction SkP. We only learn from this calcnlation that the existence of bnbble solntions 
for a given nnclens depends sensitively on the effective interactions we nse. This is to be 
expected given the fact that the mass and charge of even the lightest bnbble nnclei exceeds 
by far those of presently known nnclei. 

The resnlts displayed in Figs. 22-25 are clearly commented in the hgnre captions. 
We note in addition that the s.p. potential V{r) obtained (see Fig. 23) is intermediate 
between a sqnare well and a harmonic oscillator. Fnrthermore, we see that the energy 
density h{r) shown in Fig. 25 drops more smoothly to zero at the inner and onter snrface 
than it shonld for the leptodermons expansion to be valid. This demonstrates that the 
LDM shonld not be expected to yield a very precise valne of the bnbble energy, given the 
fact that the diameter d of the bnbble layer is in general small. 

Beside the absolnte minimnm of the energy as a fnnction of the hole parameter /, 
one may hnd one or two competing eqnilibrinm solntions at higher energy^. In Fig. 26 
we show resnlts on the system with Z=288 and iV=462. The lowest solntion is obtained 


rd 


one 


for /=0.26, the 2"'“ one with an excitation energy of 5.9 MeV at /=0.175, and the 3 
with an excitation energy of 25 MeV at /=0.12 . The difference between the total energy 
E and the LD-energy in Fig. 26 represents the shell energy. It is seen that the gronnd 
state solntion is prodnced by a strongly negative shell energy in the region of / >0.24. 
Two other solntions are mainly prodnced by the appearance of a positive shell energy 
in the region of 0.13> / >0.16 . By inspection of the level schemes (see esp. Fig. 7) 
one realizes that the strongly negative shell effect for / >0.24 is connected with a zone 
of rednced level density for Z=288 below the bnnch of rising levels with 1 radial node. 
On the other hand, in the region 0.13> / >0.16 the proton nnmber 288 and the nentron 
nnmber 462 correspond to Fermi energies which are located within the bnnch of rising 
levels with 1 radial node. 


4 Summary and discussion 

We investigated the shell correction energy of spherical bnbble nnclei. We fonnd strong 
shell effects which may give rise to shell energies of np to -40 MeV for certain magic 
nnmbers. By calcnlating the LD-energy for deformed babbles and assnming that the shell 

^We owe this insight to J.F. Berger and J. Decharge who found in their HFB calculations that up to 3 
different bubble solutions may exist for the same number of neutrons and protons. For the solution with 
the smallest value of /, these authors find the hole to be occupied by a reduced nucleon density. 
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effect due to the spherical symmetry disappears as a function of the deformation in about 
the same way as for normal compact spherical nuclei, one hnds that the fission barriers 
are of the same order of magnitude as the shell-correction energy for the spherical bubble 
solution. In favorable cases of magic numbers, this is sufficient to reduce the probability 
for spontaneous hssion practically to zero. 

The investigation shows that promising candidates for bubble nuclei are found for 
mass numbers A > 450, with the proton density being the larger the lighter the nucleus. 

The origin of the shell effects is the high degeneracy of orbitals with large angular 
momentum on the one hand, and the rapid energetic increase of radial oscillation modes 
as a function of increasing bubble radius, on the other. The details of the results (for 
instance the lower limit of mass nnmbers, where bubble nuclei begin to exist or the precise 
valne of magic numbers) depend sensitively on the choice of LD parameters and on the 
choice of the effective interaction, resp. 

The spin-orbit splitting for bubble nuclei is smaller than for ordinary nuclei and 
the sign of the splitting may be reversed. The reason is that the inner and onter surface 
regions of the nnclear density distribntion contribnte to the splitting with opposite signs 
(see Fig. 24). 

/5-decays will tend to drive the system along isobaric lines to the composition with 
the largest total binding energy. For a nnmber of systems this will lead to desintegration 
by hssion after a series of /5-decays. For systems where the shell effects vary strongly 
as a fnnction of N and Z for given A, the /5-decays will end up in a /5-stable bnbble 
composition. 

Furthermore, a bubble nucleus consisting of N nentrons and Z protons may decay 
by emission of a neutron (g=l), a proton (g=2), or an a-particle {g=3), whenever these 
decay channels are open, i.e. whenever the corresponding Q^-values are positive. Semi- 
classically, the decay probability per time nnit Wq{Qq, l;N, Z) is given by the prodnct of 
the penetrability Pq{Qq, 1) of the potential barrier 17q(r, 1) and a qnantity which can 
be classically interpreted as the number of collisions per time nnit of the particle g with 
the potential wall; 

Wq=UqPq, (4.1) 

P,(g„/) , (4.2) 

r2^ I2M 

j, = {Uq{r; /, Zq) - Qq) . (4.3) 

Here, Mi( 2 )= mass of the nentron (proton) M, M^= mass of the a-particle, Qq are the 
g-valnes of the decay reactions, and Uq represent the potential barrier between the inner 
(r|^^) and onter tnrning points (see Fig. 27): 


u. 


Vq{r) 


2Mqr^ 


ZqjZ - Zq)el 

r 


(4.4) 


Here, Vq{r) is the average nuclear potential felt by the particle g, Zq is its charge, Z is 
the charge of the mother nucleus, and I is the orbital angular momentum of the emitted 
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particle. The turning points are the solutions of the equation 

U,(r'-A) = Q, . 


(4.5) 


The frequency ujq of assaults of the barrier may be estimated from the average velocity of 
particle q in the potential well prior to emission. 

The potential V 3 (r) which acts on the composite a-particle is a rough phenomeno¬ 
logical approximation of a non-local and energy-dependent potential. We calculated the 
penetrability for an a-particle through the barrier of and the barrier of 

the bubble nucleus consisting of 470 nucleons and 182 protons (2 times For ^ 35 ^ 

and we chose experimental Q-values (Qq,(^^® 17)=4.647 MeV and Qa(^^^Fm)=8.126 

MeV) and Saxon-Woods type potential Vq{r) with parameters adjusted by Huizenga and 
Igo ig. For the bubble nucleus ^'^°182 we used the Q-value 17.7 MeV which was obtained 
from our TF-calculation based on the Skyrme-functional and the special Skyrme inter¬ 
action SkP (see Figs. 22 - 25). In all the cases the angular momentum I was chosen to be 
0. The barriers and Q-values for 3 cases of a-decay are shown in Fig. 27. We obtained 
the following values of the a-decay half-life 


ln2 

^ ■ 


(4.6) 


Using the experimentally known values of Tq, = 7-10®?/ for ^ 35 ^ _ 25 4 /^ for 

we have obtained Tq, = 168s for ^^°182 by means of linear extrapolation in the (J, log{Ta)) 

plane. 

The system ^'^°182 is a relatively light bubble nucleus. Note that we obtained this 
bubble solution only for the special Skyrme interaction SkP. As the mass and charge 
numbers of the bubble increases [A > 750, Z > 250), the Qa-values tend to become 
larger, because the binding energy per nucleon in the bubble decreases. On the other 
hand, the height of the Coulomb barrier increases and the magicity of the mother nucleus 
may reduce the Qa-value. The outcome is that the a-decay lifetime of heavy bubble 
nuclei may become much longer than the one we obtained for ^^°182. 

The bubble nuclei we studied were usually stable with respect to emission of a 
neutron or a proton, as one can infer from the negative slope ■ 

In the case of positive Q-values for p or n emission, the lifetimes can be estimated from 
the formulae (4.1)-(4.3). The decay probability Wq is smaller, the larger the angular 
momentum I of the emitted particle. Since rather large values of I occur in a typical 
bubble nucleus, the additional hindrance of the decay by the centrifugal barrier may 
become considerable. 

Another important question is to investigate the stability of a given bubble nucleus 
as a function of temperature T and the angular momentum I of the system. 

We can make the following qualitative remarks: As the temperature increases, the 
shell effect, which is essential for protecting the spherical bubble from hssion decreases 
and is known to vanish more or less completely for T > 1.5 MeV. As far as LD part of the 
energy is concerned the main effect of a hnite temperature is reduce the surface tension. 
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This increases the effective fissility parameter Xq (see Eq. (3.8)) and thus tends to favour 
the spherical bubble solution as compared to the energy of the compact spherical drop 
(see Fig. 2). On the other hand, the reduced stiffness of the surface implies that the LD 
energy of the bubble decreases faster as a function of deformation (see Fig. 21). As a 
result, bubble nuclei cannot be expected to sustain a lot of excitation energy. 

Concerning the dependence on the angular momentum again two effects have to be 
distinguished: The shell effect becomes a function of the rotational frequency Urot of the 
bubble. As a consequence, the maximal shell effect is expected to occur at a different 
value of the hole parameter /. Furthermore, the approximate effect of the rotation is 
to increase the total energy by the rotational energy of a rigidly rotating bubble. This 
rotational energy decreases as a function of /. Consequently, we may expect that an 
increase of the angular momentum of rotation lowers the energy of bubble valleys as 
compared to the energy of a rigidly rotating compact sphere. In specific cases, this effect 
may even produce a bubble valley where there would be none for angular momentum /=0 
in very much the same way as this is the case for the 2”^ valley in the region of rare earth 
nuclei. However, one should realize that the rotational energies for bubble nuclei are the 
smaller the larger the mass number A of the bubble nucleus. Consequently, for heavy 
long-lived bubbles we may not expect a large effect of the rotation. Another implication 
of the rotation will be that the optimal shape of the bubble will be deformed rather than 
spherical. 

We have not yet taken into account pairing in bubble nuclei. For nuclei away from 
shell-closures we expect the pairing energy to be generally larger than for ordinary nuclei, 
because the density of levels in the vicinity of the Fermi energy is larger due to the high 
degeneracy of single particle states of large angular momentum. The pairing energy 
matters only if at least one sort of nucleons forms a closed shell, because otherwise bubble 
configurations are unstable. This means that the pairing in the bubble nucleus is expected 
to play a role only for the non-magic sort of nucleons. We do not expect that the pairing 
leads to a qualitatively important change of the stability. 

It is of particular interest to investigate carefully the lower limits of mass- and 
charge number where sufficiently long-lived bubble nuclei can be expected. In view of 
the uncertainties concerning our knowledge of effective interactions, these predictions will 
always remain somewhat uncertain. Nevertheless, a lot of useful work can still be done. 
Assuming that one tries to produce a bubble nucleus of Nh neutrons and Zi, protons by 
a complete or incomplete fusion reaction 

(TVi, Zi) + (7V2, Z 2 ) = {N,, Z,y + ■■■ (4.7) 

the following aspects are important: 

(i) As the compact spherical nucleus with about the same neutron and proton numbers 
Nb, Zb is violently unstable with respect to fission and other decay channels, one 
may expect only a tiny fraction of the process to lead to an excited bubble nucleus. 

(ii) Assuming that AT is the difference between the initial relative kinetic energy in 

O O 

the entrance channel and the energy ZiZ 2 e^/{Ri + R 2 ) of the Coulomb barrier 
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o 

{Ri= radius of the incident nucleus (iV,, Zi)), the total excitation energy E* of the 
bubble nucleus can be estimated from the difference of the average binding energy 
per particle B{Ni,Zi) (< 0) of incident nuclei as compared to the average binding 
energy per particle B{Nb, Zb) of the bubble nucleus 

E* = [Ai ■ 5(iVi, Zi) + ^2 ■ B{N2, Z 2 ) - A ■ B{Nb, Zb)] -AE + AT (4.8) 


Here AE is the energy carried away by the other hnal particles in the reaction (4.7). 
For typical bubble nuclei, the binding energy per particle is about -4 MeV as com¬ 
pared to -7 MeV for a heavy actinide nucleus. The bracket [ ] on the r.h.s. of (4.8) 
is then in general a negative number which is to be counterbalanced by the choice 
of the kinetic energy AT. Thus the energy of the beam particles must be chosen 
considerably higher than the Coulomb barrier in order to overcome the difference 
in binding energy of the constituants of the bubble as compared to ordinary nuclei. 


We draw attention to the fact that for all examples of bubble nuclei we studied, the 
charge is “overcritical” [^. This means that the binding energy of the lowest electronic 
state is larger than 2meC^. Consequently, whenever there is a K-shell vacancy in the 
nascent bubble, it will be spontaneously hlled by an electron with simultaneous emission 
of a positon. 
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5 Appendix A1 

Eigenfunctions of the shifted square well 


The eigenfunctions (pni corresponding to the potential (2.11) have the form 


^nl 


fiiM 

andiiar) + hniyi{ar) 


for 0 < r < i ?2 
for R 2 <r < Ri 


{All) 


where /z(/3r) represents the spherical Bessel function ji{l3r) or the modihed spherical 
Bessel function ji{j3r) depending on the sgn{eni +Vi)-. 


fiiPr) = 




if ^nl +V 2 > 0 
if Enl +V2<0 


{A1.2) 


JliPr) := 

jiiPr) := ^^Ii+i{i3r) 

The functions J;_,_i(/3r) and i(/3r) are the Bessel functions and the modihed Bessel 
functions, respectively (see Ref. |]^, Chap. 9 and 10). The parameter [3 is dehned by 

1/2 


/9 = 


2M| +R 2 I 


(/ 11 . 3 ) 


The dimensionless amplitudes Uni, hni and the eigenenergies e^i are determined by the 
boundary and continuity requirements 


fi{(3R2) = aniji{aR2) + bniyi{aR2) 
(3f{{f3R2) = a [anij'i{aR 2 ) + bniy'i{aR 2 )] 
0 = ttnijiiaRi) + bni yi{aRi). 

Here and in what follows we use the notation: 

'dji{z)\ 


(H1.4) 

(H1.5) 

(H1.6) 


fi{aR2) = 


dz 


and equally for y[{aR 2 ), // {PR 2 ) ■ 


z=aR2 


We satisfy the eqns. (Al.6) by the ansatz 

(^nl (^nl yi{c^Rl) 
bill 

The amplitude is then obtained from (Al.4) 

fi{f3R2) 


0>nl — 


yi{aRi)ji{aR2) - ji{aRi)yi{aR2) ' 
and the eigenenergies Sni from the continuity of the logarithmic derivatives 

yi{aRi)j[{aR2) - ji{aRi)y[{aR2) 

= Cl 


(Al.l) 

(Al.S) 

(AIJ) 


MPRi) yi{aRi)ji[aRi) - ji(aRi)yi(aRi) ' 

The overall normalization constant Nni is of course dehned by the condition (2.10). 


(Al.lO) 
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6 Appendix A2 

Eigenfunctions of the shifted harmonic oscillator 


For r < R 2 , the radial wave functions are again given by (Al.l) and (Al.2) up to a 
factor Cni (see (A2.4)). At the point r = i? 2 , the radial wave function (pni(r) in region I 
((phi(r)) and in region II must be continuous and also its derivatives: 

dR2 dR2 

Furthermore, y:>hi{r) must satisfy (2.5.2) at infinity. 

The eigenvalues e^i and the wave functions ^hi{r) must be determined numerically 
or by the WKB approximation. In both cases it is convenient to introduce the radial 
functions ul^i{r]) depending on the dimensionless variable r] 


(212.1) 

{A2.2) 


Tj = jr = 

for describing the functions (pniif) in region I: 


'Mo; 


h 


(A2.3) 


, , . CnifiiPr) for 0 < r < 

^ni{r)-Nni{ for r>i ?2 . 


{A2A) 


Henceforth, we shall suppress the subscripts nl for simplicity. The radial function u^{ri) 
satishes the differential equation 


d‘^u^{r]) 

drf 


l{l + l) 

rj2 


where 


K : = 


+ {r]-r]Y + K 

fj = ■jR 
2{e +Vo) 


u\v) = 0, 


huj 


(A2.5) 

(A2.6) 

(A2.7) 


For large values of rj, the centrifugal term in (A2.5) is negligible and the equa¬ 

tion (A2.5) assumes the standard form of the differential equation for parabolic cylinder 

(A2.8) 


functions (see Ref. [j^, p. 686): 

d?w 




when introducing the new independent variable. 


^ := V2{7] - 7]). 


(A2.9) 
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(A2.10) 


The functions v}{rj) approach the solutions w{^) of Eq. (A2.8) 

u^iri) i) = w[k] V2{t] - fj)] 

for r^-values satisfying the inequality 

77^(77 - >/(/ +1). (A2.11) 


For obtaining normalizable eigenfunctions (pni{f') we need solutions w{^) of the differential 
equation (A2.8) which go to zero at inhnity (.^ —*• cxo) at least ~ The mathemat¬ 
ical properties of the two basic, linearly independent solutions 7 / 1,2 of Eq. (A2.8) 

Eq. (A2.13) and (A2.14) are described in Chapter 19 of Ref. 
combination which vanishes for > cx) is given by 


10 


see 

The special linear 


W(K, 


;0 = 0F24 


yiiO 247/2(0 


2ir(^) r(i^)_ 


7/i(0 = e-wM 


1/2(0 = ee--M 


4 ’2’ 4 J 

'3-k 3 0\ 

^^r’2’Tj ■ 

Here, M(a, b; z) is the confluent hypergeometric function 

MK6;.) = l + -- + ^- + 


(A2.12) 

(A2.13) 

(A2.14) 


(A2.15) 


It is interesting to consider the special case of k = 277 -|- I (where n = 0,1,2,...), where 
the functions w{k]^) become equal to the even Hermite functions for even n and to the 
odd Hermite functions for odd n (see Ref. |^, Eq. (19.13.1)) 


w{k = 2n + 1]^) = e 4 ^^ He^ 


(0 


(A2.16) 


,2 ( p2 'I 

//e00:=(-ire^ —|e-^|. (A2.17) 

i.e. to the eigenfunctions of the linear harmonic oscillator. Of course, the functions 77 ^( 7 /) 
do not even asymptotically tend to this limit because of the boundary conditions (A2.1-2) 
at r = R 2 and because of the centrifugal term in (A2.5). In the case of a general value 
of K, the functions w{k-,^) represent ’’parabolic cylinder functions” and are frequently 
denoted by (Ref. [jT^, Eq. (19.3.1) 


=Pi(^_i)(0- 
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(A2.18) 

















For 1^1 3> -Y) they have the following asymptotic form (Ref. 0, (19.8.1)): 


^ {k — 1)(k — 3) ^ {k — 1)(k — 3){k — 5){k — 7) 


. 23^2 25 - 4^2 

(712.19) 

Solutions u^{ri) which tend asymptotically to the parabolic cylinder function (see 

(A2.10)) are thus normalizable for all values of k. 

The precise values of k (or Sni) and of the amplitude ani are determined by the 
boundary conditions (A2.1-2). Substituting (A2.4) these boundary conditions read 


Cnlfim2) = 


Mnz(72^2) 


7 -R 2 


CnZ/3 


'dfijz) 

dz 


I 


= 7 ■ 


Z=PR2 


dr} 




V 




(A2.20) 

(A2.21) 


r)=7i?2 


The functions u\i{ri) can be determined numerically by the ’’shooting method”. 


7 Appendix A3 

Thomas-Fermi approximation for eigenenergies £ni and single¬ 
particle densities Pni{r). 


We mentioned already the WKB-approximation for the unperturbed eigenvalues Sni 
in Chapter 2. The Eq. (2.15) represents the simplest version of the WKB-approximation 
which contains no correction from the vicinity of the turning points nor from the classically 
forbidden regions. In what follows we only refer to this simplest version and claim that 
it can essentially also be obtained from the Thomas-Fermi approximation (TFA). 

Let us hrst assume that a given eigenvalue Sni has been calculated from the simple 
semiclassical quantization condition (2.15) which we now write in the form 





hn, 


where the potential Ui{r) is dehned by 


yir) + neutrons 

V(r) + eoVcb(r) + for protons 


(A3.1) 


(A3.2) 


and where n = 0,1, 2 ,.... 

As we stated already in Sect. (2.1), the Coulomb potential Vcb(A) produced by the 
density distribution Pp(rJ of the protons (see Eq. (2.16)) can be easily evaluated in closed 
form if we approximate the density Pp(r) by a uniform distribution in the bubble layer 
(see (2.18)). The result is given in Eq. (2.19). 
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The spin-orbit potential Vso ((2.33’) or (2.37)) acting on eigenstates (2.2) in the 
coupled representation is a c-number operator. Consequently, in order to evaluate its 
mean value, we only need the single particle densities p^i (r) as defined in Eq. (2.40). 
Since no derivatives occur, the accuracy requirements concerning these single-particle 
densities are less stringent. Therefore, we determined these densities in the Thomas- 
Fermi approximation (TEA). 

The numbers N of neutrons and Z of protons are related to the Fermi energies 
and Ep by the Eqs. (2.25). On the other hand, the total densities (r), Pp (r) of n,p 
must satisfy the trivial relations 

N = jdh p, (r) , (A3.3) 

Z = J d\ Pp (r) . (ASA) 

By comparison with (2.25) we obtain the TEA of the densities: 





(A3.5) 


Pp [r) = 




d^pOo\e^F-^-V{r)-Vc^{r) 


87r(2M)3/2 

3/l3 


[e^F-V{r)-Vcy>{r)f\ 


(A3.6) 

By integrating over the phase space enclosed by two neighboring surfaces corresponding to 
the constant energies (^e^, — and (ei, + we obtain the number of n or p contained in 
this phase space. Requiring that the volume is occupied by 1 n or p we may determine the 
corresponding energy width Ae which depends on the energy Ei, and may be insignificantly 
different for n and p: 


1 = 


Ae 


p 




d\- d^pOo + 


+ + (A3.7) 


2 2M 

\ / \ / 

Here V{r) is the potential acting on a n or p, resp.: 

V{r) = V{r) for neutrons (A3.8) 

V{r) =V{r) + eoVcb(^") for protons. (A3.9) 

The density distribution p (r; Ei,) of the neutrons (protons) in the energy layer is given by 


P (r;e^) = ^ f d^pdo (e^ + ^ 


or 




f, (r;E„) = 


P 

2M 


V{r) 9, 


p. 


A5' 


+ - . (A3.10) 


Ey -\- 


Ae 


V{r] 


3/2 


Ae 


V(r] 


3 / 2 ' 


(A3.11) 
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where the width Ae of the energy layer is to be determined from relation (A3.7). The 
condition (A3.7) reads more explicitly 


1 




3/2 




(A3.12) 


In each case, integrations over the radial coordinate r extend over the interval where the 
integrand is real. 

For calculating the average of the spin-orbit potential l4o we need the single particle 
density (r) of a nucleon of given energy Sni and given orbital angular momentum /. 
The simple form (A3.11) of the density cannot be used because it averages over particles of 
all angular momenta /. From the TFA we cannot obtain a single particle density with the 
desired number n of radial nodes but we may obtain an approximate density distribution 
for a nucleon of given orbital angular momentum I by performing the integration in (A3.7) 
and (A3. 10 ) under the constraint 

£2 = (r X:p)2 = /(/ + l)h^ (A3.12) 

This can be simply achieved by choosing the 3-axis of the momentum integration in the 
direction of the position vector r and using cylinder coordinates in p-space: 


P = P3er+p ±(cos Cl - 1 - sin 62 ), 
d^p = dpz dp±p±d(j)p , 

r 2 r 2 2 2 

L = L = r pj_ 

r2 

P = P3+P±= P3 + ^ ■ 

The condition (A3. 12 ) can be fulhlled conveniently by introducing the continuous dimen¬ 
sionless integration variable A instead of p± 

L Xh 


and inserting the (5-function 6 \^X — + l)j into the integrals. We thus obtain the 

following 2 relations in analogy to the Eqs. (A3.7) and (A3.10) 

1 = %Sdh-Sdp,dXdct>p^,e, M _ |L _ ^ _ ^(p)) 

■ ^0 (tw + - ^/(/ + 1)) 

Pnl (^; ^nl) = fdps dX d<i)p X 00 (enl + ^ ^ ^ W ) .^3 

■ ^0 + v{r) - eni + (5 (a - + 1)) . 
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Carrying out the integrations in (A3.14) we obtain 


Pnl ir 1 ^nl) 


2V2M. 


2Tihr‘^ 




(A3.15) 


where Ui{r) represents the potential including the centrifugal term (see Eq. (A3.2)). For 
/ 3> 1, we have 2^l{l + 1) 21. This represents the classical approximation for the 

quantum-mechanical degeneracy factor ((2/ -|- 1). The relation (A3.13) takes the form: 


2Jl{l + 


drlJeni + ^- Ui{r) - ^ “ ^lir-) 


(A3.16) 


The main dehciency of the TEA is to completely neglect the density in the classically 
forbidden region and to ignore the correct nodal structure of the density distribution of 
single particles inside the potential well. 

The TEA is closely related to the WKB approximation. Let us consider the Hamil¬ 
tonian 

(A3.17) 

of a 1-dimensional potential model and the corresponding classical Hamiltonian function 




(A3.17') 


If the particles moving in the potential V{x) are fermions of spin the number N{eu) of 
fermions which £11 the energy levels below a given energy are given in the TEA by 


= J dp9o{e^-H{x,p)) . 


(A3.18) 


We can use the relation (A3.18) in order to determine approximate eigenvalues ei < 62 < 
... by requiring M^ey) to be equal to 2z/, where v = 1,2 ,... counts the eigenvalues in 
rising order. We thus obtain the implicit equation for determining the eigenvalues 


dx J dp9o (sy — H{x,p)) . 


(A3.19) 


For the linear harmonic oscillator V{x) = we find from (A3.19) the eigenvalues 


Ey = VhUJ 


(A3.20) 


for V = 1,2,.... Using the same definition of the index z/, the correct eigenvalues are 
given by 

ey = {i'-^hu. (A3.21) 
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For the 1-dimensional infinite square well, the TEA (A3.19) yields the correct result. In 
the case of a spherical potential in 3-dimensional space we have to restrict the integration 
of the phase space to the hypersurface defined by the constraint (A3.12), i.e. we apply 
the method for calculation of the density distribution This yields 

Ti^ ■ + 1 ) f ^ f ( pi \ 1 

n-i 2 l + l) = -- J dr J dp^e ^ - Uiir)j ^ , (A3.22) 

where n = 1,2,... and {21 + 1) is the degeneracy of a level of given /. Within our 
approximation we should put 2 ^Jl{l + 1) ^ 2/ ^ 2/ -|- 1. Consequently, (A3.22) takes the 
form ^ 

V 2 M f drJeni — Ui{r) = nTm , (A3.23) 

J a 

where a and b are the lower and upper turning points and n = 1, 2,.... The Eq. (A3.23) 
resembles the WKB formula given by Eq. (2.15). We note, however, that in (A3.23) the 
radial quantum number n starts from 1 whereas in (2.15) it starts from 0. For n ^ 1, the 
results are equivalent. 
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Figures captions. 

1. The function V 2 {Xq) defining the stationary points of the energy S{v 2 ]Xq) as a 
function of the fissility parameter Xq. 

2. Difference S = XEu^/Es{Ro) of the LD energy of a bubble of reduced inner radius 
V 2 and the LD energy of a compact spherical nucleus of radius Rq. 

3. Lines of equal bubble fissility parameter X in the (iV, Z) plane. Isobars are indicated 
by the parallel straight lines. 

4. Lines of eqnal binding energy gain AEi^d by bnbble formation. Parallel straight 
lines indicate isobars. 

5. Lines of eqnal hole fraction / in the {N, Z) plane. Parallel straight lines indicate 
isobars. 

6 . Lines of constant LD-binding energy per particle. Parallel straight lines indicate 
isobars. 

7. Level scheme for the shifted inhnite sqnare well with the spin-orbit conpling as a 
fnnction of the hole fraction / = {R 2 /R 1 Y (see Eq. (3.15)). The letters s, p, d, 
f, g, h, i, j, ... mean the orbital angnlar momenta I = 0,1,... . For 3 valnes of 
/ (0.12, 0.24, 0.28) nnmbers are written jnst above certain single particle energies 
El, whenever the distance {e^^i — e^) exceeds 1.5 energy nnits. They represent the 
nnmber of nentrons (or of protons) occnpying single particle states with energy 
Sk < Su- 

8 . Level scheme for the shifted harmonic oscillator with spin-orbit conpling as a fnnc¬ 
tion of the hole fraction /. The oscillator freqnency uj was determined as a fnnction 
of / for A=500. 

9. Level scheme for the inhnite square well potential with a hnite step in the inner 
region. Height (Vq ~ ^ 2 ) of the step: 60 MeV for r < i? 2 - 

10. Same as in Fig. 9, bnt with a step height of 40 MeV. 

11. Shell correction energy SEgheii for bnbbles with /=0.24 as a fnnction of iV or Z for 
the shifted inhnite sqnare well with spin-orbit term. 

12. Same as in Fig. 11, bnt for a sqnare well with an internal step height (Vq — 7^) = 
40, 50, and 60 MeV. The case of the simple inhnite sqnare well is given again for 
comparison. 

13. Same as in Fig. 11, bnt for a shifted harmonic oscillator potential. The oscillator 
freqnency corresponds to the mass nnmber H=500. 
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14. Shell energy evaluated from the inhnite square well plus I ■ s'-term as a function of 
N (or Z) and the hole fraction /. The energy units are [100/A^/^] MeV. 

15. As in Fig. 14, but for a shifted harmonic oscillator potential corresponding to the 
particle number y4=500. Energy unit: MeV. 

16. As in Fig. 14, but for a shifted harmonic oscillator potential corresponding to the 
particle number A=1000. Energy unit: MeV. 

17. Lines of constant shell energy as a function of particle number [N or Z) and the 
hole fraction parameter / for a shifted inhnite square plus I ■ s-coupling. The energy 
unit is [100/A^/^] MeV. 

18. As in Fig. 17, but for a shifted harmonic oscillator with / -s-coupling. The oscillator 
frequency corresponds to the particle number A=500. Energy unit: MeV. 

19. As in Fig. 18, but with the oscillator frequency chosen for the mass number A=1000 

20. The map of the liquid drop plus the shell correction energy corresponding to the 
bubble solution. 

21. LD-energy (solid line), shell correction energy (dotted line), and the total energy 
(dashed dotted line) as a function of the quadrupole deformation (52 of the outer 
bubble surface Si. The LD-energy is minimized with respect to the multipole 
deformation of the order 4 and 6 of the surface Si and the multipole deformation 2 
to 6 of the inner surface 5*2. The long dashed line represents the LD-energy when 
both surfaces have the same deformation. The short dashed line corresponds to the 
case when the inner surface is a sphere. 

22. Density Pp{r) of protons (solid line) and Pnij") of neutrons (dashed line) obtained by 
variation of the Skyrme functional (3.28) obtained within the EFT approximation 
for SkP forces and for the nucleus Z=184, A=470. 

23. Nuclear central s.p. potential acting on protons (solid line) and neutrons (dashed 
line) as obtained from the variation of the Skyrme energy functional H- 

24. Average spin-orbit potential Visi^) A=470 and Z=184. 

25. Energy density h{r) for A=470 and Z=184. 

26. Total energy E and liquid drop energy 77^^ as a function of the hole parameter / for 
a system with Z=288 protons and A=750 nucleons. 3 bubble solutions correspond¬ 
ing to the /-values /=0.26 (ground state), /=0.175 (bubble valley at the excitation 
energy E*=5.9 MeV), and /=0.12 (bubble valley at the excitation energy E*=25 
MeV) are seen to exist. The inner and outer radii of the 3 solutions are: i?2=7.69 fm 
and i?i=12.05 fm for /=0.26, i?2=6.50 fm and i?i=11.62 fm for /=0.175, i?2=5.61 
fm and i?i=11.37 fm for /=0.12 . 
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27. The potential barrier 7/3 (r) acting on an a-particle in ^^^77, in and in the 

bubble nucleus ^'^°184 are shown together with the corresponding Qs = Qa values. 
The bubble nucleus ^^'^184 is the one obtained as a result of a TF-calculation with 
the Skyrme interaction SkP (see Figs. 22 - 25). 
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